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Abstract 

A experimental study of low Reynolds numbers unsteady pipe flows of a yield stress shear thinning fluid (Carbopol® 
- 980) is presented. The investigation of the solid-fluid transition in a rheometric flow in the presence and in the ab- 
sence of the wall slip reveals a coupling between the irreversible deformation states and the wall slip phenomenon. 
Particularly, the presence of wall slip nearly suppresses the scaling of the deformation power deficit associated to the 
rheological hysteresis with the rate at which the material is forced. The irreversible solid-fluid transition and the wall 
slip behaviour emerge in the same range of the applied stresses and thus, the two phenomena appear to be coupled 
to each other. In-situ measurements of the flow fields performed during an increasing/decreasing stepped pressure 
ramp reveal three distinct flow regimes: solid (pluglike), solid-fluid and fluid. The time resolved measurements of low 
Reynolds number unsteady flows driven by an increasing/decreasing controlled pressure ramps reveal a irreversible 
yielding scenario similar to that observed in the rheometric flows. The deformation power deficit associated with the 
hysteresis observed during the increasing/decreasing branches of the pressure ramps reveals a dependence on the rate 
at which the unsteady flow is driven consistent with that observed during the rheological measurements in the presence 
of slip. The dependence of the slip velocity on the wall shear stresses reveals a Navier-type slip behaviour only within 
the fluid flow regime, which indicates that the wall slip phenomenon is directly coupled to the solid-fluid transition. 
A universal scaling of the slip velocity with the wall velocity gradients is found and the slip length is independent on 
the characteristic time of forcing to. The paper closes with a discussion of the main findings, their possible impact on 
our current understanding of the yielding and slip behaviour of Carbopol® gels. Several steps worth being pursued 
by future experimental/theoretical studies are proposed. 
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fo = 15s, (b) fo = 20s. The squares (□) refer to the increasing branch of the stress ramp and the circles 

(o) refer to the decreasing branch of the stress ramp |l5] 

7 Dependence of the area of the hysteresis of the dependencies of the mean flow velocity on the driving 
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eye,Aocf- 01 |l6| 

8 Dependence of the slip velocity U s on the applied pressure drop Ap for various values of the charac- 
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Dependence of the slip velocity U s on the wall velocity gradients. The full line is a guide for the eye, 

s A) - 1 = 4 (j), (o, .) - 1 = 7.5 (j), (□, ■) - r = 10 (j), 
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(V, A) - 1 = 15 (s), (<]. O - ? = 20 (i) 

10 Dependence of the slip velocity t/ s on the wall shear stresses T w . The full line is a power law fit, U s = 
2.3 ■ 10~ 4 (±5.4 ■ 10~ 5 ) ■ T v 5/ 32±0 ' 26 . The dotted and the dash-dotted lines are the theoretical predictions 
by Piau (see Ref. |Piau| fl2007[ )): U s °c t v 1/ 3 , U s °< t 2 . The symbols are: (A, A) - 1 = 4 (s), (o, •) - 
fo = 7.5 (j), (□, ■) - f = 10 (s), (V, A) - f = 15 (s), (<, +)-to = 20 (s) 
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1. Introduction 



Yield stress fluids represent a broad class of materials made of high molecular weight molecular constituents 
which display a solid-like behaviour as long as the stress applied onto them does not exceed a critical value called the 
yield stress, T v , and a fluid one beyond this threshold. 

The constantly increasing level of interest of both theoreticians and experimentalists in yield stress fluids has, in 
our opinion, a two-fold motivation. From a practical perspective, such materials have found an increasing number 
of practical applications for several major industries (which include foods, cosmetics, oil field, etc.) and they are 
encountered in the daily life in various forms such as food pastes, hair gels and emulsions, cement, mud and more. 
Carbopol® gels have been widely used in many rheological and hydrodynamic experimental studies and they have 
been considered for over a decade a model yield stress fluids. The wide use of Carbopol® gels is equally due to 
their optical transparency which allows an in-situ assessment of the flow fields, their stability (they can be stored 
for extended periods of time without noticeable changes in their rheological behaviour), [Piau| (2007). Carbopol® 
resins are synthetic polymers of acrylic acid initially introduced over six decades ago (B.F. Goodrich Co.). They 
are cross-linked with various chemical compounds such as di vinyl -glycol, allyl-sucrose, polyalkenyl polyether etc. 
In an anhydrous form, the average size of the polymer molecules is of the order of hundreds of nanometers. In the 
absence of cross-links, the polymer particle can be viewed as a collection of linear chains intertwined (in a coiled state) 
but not chemically bonded. The Carbopol® micro-gel particles are soluble in polar solvents and, upon dissolution, 
each individual polymer molecule hydrates, partially uncoils, and swells several orders of magnitude. Addition of 
a neutralising agent leads to the creation of negative charges along the polymer backbone due to ionisation of the 
carboxylic acid groups. Consequently, swollen polymer molecules cross-link, forming a system of micro-gel particles . 
The micro-gel system can sustain finite deformations (behaving like an elastic solid) prior to damage. When local 
applied stresses exceed a threshold value the gel system breaks apart and the material starts to flow. This is the 
commonly accepted microscopic scale origin of the macroscopic yielding of the material. 

From a fundamental perspective and even in the case of minimal thixotropic effects, yield stress materials con- 
tinue triggering intensive debates and posing difficult challenges to both theoreticians and experimentalists. From a 
historical perspective, the best known discussion is related to the very existence of the yield stress Barnes ( 1999| ): 
Barnes and Walters (1985). It was thus suggested that the yield stress is a engineering reality rather than a physical 



one or, in other words, the apparent yield stress behaviour observed during rheological tests is related to a very large 
Newtonian viscosity rather than to the onset of flow (yielding). 

In context of the recent developments of the rheometric equipment (now able to resolve torques as small as 
Q.lnNm and rates of shear as small as 10~ 7 s~ l ), a further pursue of this already classical debate on the existence of 
a yield stress is, in our opinion, trite. Indeed, recent rheological tests performed on various Carbopol® gels proved 
unequivocally the existence of a true yielding behaviour, Putz and Burghelea ( 2009| ); [Vloller et al.| ( |2009[ ); Bonn and 
Denn ( |2009| ); penn and Bonn] ( |2011[ ). 

The existence of a Hookean solid deformation regime preceding the solid-fluid transition has been experimentally 



demonstrated in Ref. Putz and Burghelea (2009) in the form of a plateau of the rate of deformation measured during 



increasing/decreasing linear stress ramps. The following characteristics of the transition from a solid like to a fluid like 
behaviour in a rheometric flow (a plate-plate geometry) in the absence of wall slip have been found experimentally: 

1 . The solid-fluid transition is not direct but mediated by an intermediate deformation regime characterised by 
a coexistence of the solid and fluid behaviour. This experimental facts render the classical Herschel-Bulkley 
formalism inapplicable over the entire range of material deformations. 

2. Except for the fluid deformation regime observed for applied stresses significantly larger than the yield stress 
Tj , the solid-fluid transition is not reversible upon increasing/decreasing the applied stress and a non-negligible 
hysteresis effect is observed. 

3. In the absence of wall slip, the deformation power deficit associated to the experimentally observed hysteresis 
scales as a power law with the rate at which the material is deformed, i. e. the inverse of the time interval the 
applied stress is maintained constant during the increasing/decreasing branches of the stress ramp. Therefore, 
only in the asymptotic limit of very slow deformations the hysteresis tends to disappear and the solid-fluid 
transition becomes reversible. 
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The experimentally found features of the solid-fluid transition in a Carbopol® gel enumerated above cannot be 
captured within the classical theoretical frameworks (Bingham, Herschel-Bulkley, Casson) which predict a reversible 
deformation and a well defined onset of flow. 

The rheological characterisation of several structurally and physico chemically different yield stress fluids (a 
Carbopol® gel, a commercial hair gel, a commercial shaving foam and a cosmetic water-in-oil emulsion) presented 
by M0ller an his coworkers bring additional strong arguments for the existence of the yield stress and seems to 
rule out the possibility of a large Newtonian viscous behaviour prior to yielding, M0ller et al. (2009). Whereas a 



viscosity plateau is indeed observed within a range of low applied stress, its value increases as a power law with the 
characteristic time to during each the shear rate data is averaged for each value of the applied stress, 770 °= (with 
a ps 1). Thus, the authors conclude that the apparent viscosity plateau is only due to the fact that a steady state of 
deformation is not achieved for low applied stresses and therefore it should be interpreted with caution (see Fig. 2 
in Ref. M0ller et al. (2009)). The experimentally found power law scaling of the plateau viscosity observed prior to 
yielding implies a divergence of this viscosity which M0ller and his coworkers interpret as an unequivocal proof for 
the existence of a true yield stress. 

Th e insufficiency of the c l assical descriptions of the yield stress fluids has been clearly emphasised by Denn and 



Bonn, Denn and Bonn| ( 201 l ); |Bonn and Denn| ( [2009| ). To overcome these difficulties, Denn and Bonn propose in Ref. 
|Denn and Bonn "( 2011 ) several ingredients which a realistic theoretical model should contain: thixotropy, avalanche 
behaviour, loss offore-aft symmetry inflow, shear banding without a stress gradient, shear localisation. At a technical 
level, they equally suggest employing a kinetic equation for a structural parameter which dynamica lly characterises 



Denn and Bonn 



t he am ount of un-yielded material as a function of the externally applied stress, see Eq. 7 in Ref. 
( 2011 ). To these ingredients we would also add the practical need of accounting for the wall slip behaviour. 

Part of these tasks (which in our opinion are fully justified by the compelling experimental evidence) have already 
been addressed (at least partially). A simplistic kinetic model that reveals a thixotropic behaviour (time dependent 
yielding, irreversibility of deformation states), can phenomen ologically explain the fore-aft symmetry breaking during 
flow past moving rigid body observed experimentally in Ref. Putz et al. ( 2008 ) and is co nsistent with a shear bandin g 
(phrased as coexistence between solid and fluid behaviour) has been proposed in Ref. Putz and Burghelea ( 2009| ). 
Although this model proved its ability to describe the solid-fluid transition in a rheometric flow (e.g. it can accurately 
fit stress strain data acquired during controlled stress ramps), additional theoretical developments are needed in order 
to describe more complicated flow situations such as sedimentation problems, start-up pipe flows etc. We also note that 
this oversimplified model was not formulated from first principles and it is still lacking a thermodynamic validation. 

Whereas the recent experimental evidence for the existence of yield stress seems compelling to us, we do believe, 
however, that dynamics of the yielding transition remains an interesting topic for both theoreticians and experimental- 
ists, particularly in realistic flow settings where a steady state of deformation can not be achieved and in the presence 
of slippery boundaries. In practical flow situations involving pasty materials the yielding transition is often associated 
with the wall slip phenomenon which complicates even further the theoretical description. Contrarily to the case of 
polymer melts where the slip is observed at large rates of strain Denn (2001, 2008), the wall slip in pasty materials 
emerges within a range of small strains. This implies that in the case of pasty materials the wall slip and the yielding 
transition are most probably coupled to each other which brings an additional experimental and theoretical difficulty 
in exploring the solid-fluid transition and modelling realistic flows of such materials in industrially relevant settings. 

A detailed experimental investigation of coupling between the wall slip phenomenon and yielding observed in a 



coaxial cylinder geometry for various pasty materials has been reported by Bertola and his coworkers, Bertola et al 



( P003[ ). The novelty of their study comes from combining classical rheometry techniques (increasing/decreasing con- 
trolled stress ramps) which provides the characterisation of wall slip and yielding at an integral scale with the magnetic 
resonance imaging (MRI) technique which provides local measurements of the velocity profiles. The successful com- 
bination of these complementary techniques allowed Bertola and his coworkers to provide a detailed description of 
the coupling between sli p and yielding and to identify a critical rate of shear below which no steady state flow fields 
can be observed by MRI, |Bertola et al.| ( |2003| ). 

A detailed theoretical study of the combined effects of compressibility and slip of Poiseuille flows of a Herschel- 
Bulkley fluid has been recently presented by Damianou and coworkers, pamianou et al. (2012). A central conclusion 
of this work is that, in the incompressible flow case, the flow field may become plug-like corresponding to a critical 
value of the slip parameter which is found to be proportional to the yield stress of the material. 

Whereas there exists a significant number of studies of the yielding of Carbopol® gels in rheometric flows, the 
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body of literature concerning experimental studies of the solid-fluid transition of Carbopol® gels in industrially rele- 
vant flows is, to our best knowledge, more limited. A experimental study of the yielding and flow properties of a 0.2% 
Carbopol® in a capillary tube has been performed by Gonzalez and his coworkers, |Perez-Gonzalez et~aT| ( |2012[ ). By 
combining classical rotational rheometry techniques with Digital Particle Image Velocimetry (DPIV) measurements 
of steady flows driven at various pressure drops the authors of this study conclude that a yield stress behaviour does 
exist. As the driving pressure drop is gradually increased, they observe a smooth solid-fluid transition. To overcome 
the difficulty of accurately measuring the yield stress triggered by the smoothness of the solid-fluid transition, the 
authors of the study propose an alternative method of measuring this parameter based on the analysis of the measured 
velocity profiles, Perez-Gonzalez et al. (2012). As the flows studied in Ref. Perez-Gonzalez et al. (2012) were sta- 
tionary, these experimental results can neither confirm nor contradict the main features of the solid-fluid transition 
observed during unsteady rheometric flows in Refs. Putz and Burghelea| ( |2009 ); Weber et al. (2012) which motivated 
us to consider the unsteady pipe flow case in the present study. 

The body of existing literature on unsteady pipe flows of Carbopol® gels is even more limited. Park and his 
coworkers have studied both experimentally and theoretically oscillatory flows of a 0.075% Carbopol® solution in 
a acrylic made U-shaped tube with a circul ar cross section , Park and Liij ( 201C ). Consistently with the findings 
reported in Ref. Putz and Burghele^ (2009), the authors of this study have provided a clear experimental evidence 
on the significant role of elasticity on the measured unsteady velocity profiles particularly near the centreline of the 
flow channel. Due to the elastic effects, strongly nonlinear time histories of the strain are measured. To model their 
experimental findings, Park and his coworkers propose an elasto-viscoplastic model which fits their data better than 
the regularised Bingham model. Some discrepancies between theory and experiment still exist and may be due to the 
fact that their model does not account for the velocity slip near the wall. 

The central aim of the present study is to present a detailed experimental characterisation of unsteady flows of a 
Carbopol® gel in the presence of wall slip in both a rheometric flow and a pipe flow. The general question addressed 
through our study concerns with the dynamics of the solid-fluid transition observed during the unsteady low Reynolds 
numberpipe flow of a Carbopol® gel and its connection to both the wall slip phenomenon and the solid-fluid transition 
previously observed in rheometric flows (with and without slip at the boundaries). 

The paper is organised as follows. The experimental setup, the measuring technique and the physical and Theolog- 
ical properties of the solutions are detailed in Sec. ^| The experimental results are presented in Sec. || A systematic 
description of the solid-fluid transition observed in a rheometric flow with a particular emphasis on the role of wall 
slip and its coupling to the irreversibility of the deformation states is presented in Sec. [3. 



A detailed characterisation 



of inertia-free unsteady channel flows of a Carbopol® solution is presented in Sec. 3.2 The paper closes with_a 
discussion of our main findings, of their impact on the current understanding of flows of yield stress fluids, Sec 
Several directions worth being pursued by future studies will also be discussed. 



2. Experimental setup and methods 

2.1. Experimental apparatus 

The experimental setup is schematically illustrated in Fig. [j]. It consists of a flow channel FC of length L = 1 15 cm 
with a circular cross-section of inner radius R = 1 .9 mm. 

The flow channel is housed by a acrylic made fish tank FT filled with water circulating at a constant temperature, 
T = 23"C. Besides its role in maintaining a constant operating temperature, the water circulation around the flow 
channel also provides a good matching of the refractive index which minimises the optical distortions induced by the 
curved surface of the flow channel. 

A green laser beam (A =514 nm) with a power of 500 mW emitted by the solid state laser L (from Changchun 
Industries, Model LD-WL206) is deflected by a mirror through a cylindrical optics block CO which reshapes it into 
a horizontal laser sheet, Fig. Ill The cylindrical optics block is composed of a glass rod with a short focal distance 
fi~2 mm and a cylindrical lens with a larger focal distance, m 5 cm. The two optical elements are mounted 
orthogonally to each other and in a telescopic arrangement such that the primary horizontal laser sheet generated by 
the glass rod is focused on the vertical direction in the middle of the flow channel by the second lens minimising thus 
its thickness in the measurement region. The thickness of the generated laser sheet is roughly 80/Xm in the beam waist 
region which is positioned at the centreline of the flow channel. The working fluid was seeded with an amount of 
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200 parts per million (ppm) of polyamide particles with a diameter of 60 \lm (from Dantec Dynamics). A time series 
of flow images is acquired with a 12 bit quantisation digital camera CCD (Thorlabs, model DCU224C) through a 
12X zoom lens (from La Vision). The resolution of the acquired images is 1280X1024 pixels which insures a spatial 
resolution of 2.96/ira per pixel. 

2.2. Preparation of the Carbopol® solutions 

A 0.08% (by weight) solution of Carbopol® 980 has been used as a working fluid. The procedure for the prepa- 
ration of the solution is described as follows. First, the right amount of anhydrous Carbopol® 980 has been gently 
dissolved in water while continuing stirring the mixture with a commercial magnetic stirring device. The stirring 
process was carried on until the entire amount of polymer was homogeneously dissolved. A particular attention has 
been paid to the homogeneity of the final mixture, which has been assessed visually by monitoring the refractive 
index contrast of the mixed solution. Next, the pH of the mixture (initially around 3.2, due to the dissociation of the 
polyacrylic acid in water) has been brought to a neutral value by addition of about 140 parts per million (ppm) of 
sodium hydroxide (NaOH). The final value of the pH has been carefully monitored using a digital pH- meter (from 
Grosseron). Upon the ionisation and the neutralisation of the carboxylic groups of the polyacrylic acid, the Carbopol® 
particles swell dramatically forming a percolated micro-gel, ^iaij ( 2007). As a particular attention has been paid to 



control the neutral pH, the concentration of the sodium ions was presumably equal to that of the carboxylate groups. 
To homogenize the gelled mixture, the sample has been stirred with a propeller mixer (at 400 rpm) for about two 
hours and then allowed to rest for several more hours. Prior to each rheological test, the air bubbles entrapped into the 
gel during the stirring process have been removed by placing the samples in a low pressure chamber for 30 minutes. 
Finally, the fluid batch has been stored at room temperature in a sealed container in order to minimize the evaporation 
of the solvent. 



2.3. Rheological characterisation of the Carbopol® solutions 

The rheological properties of the solutions were investigated using a Mars /// (from Thermofischer) equipped with 
a Peltier system able to control the temperature with an accuracy better than 0.1 °C and a nano-torque module which 
allows one to accurately explore very small rates of deformation, y « 10~ 5 

A first major concern for the rheological measurements was the emergence of the wall slip phenomenon at the 
contact with the measuring geometry, which has been previously recognised as a major source of experimental arte- 
facts when Carbopol® gels are used, f'iau (2007); Piau and Debiane ( |2004 ). To quantitatively assess the role of the 
wall slip, two alternative geometries have been used: a smooth parallel plate system and a serrated one. The radius 
of each geometry was R — 35 mm and the gap between the parallel plates was d = 1 mm. Measurements of flow 
curves performed with the rough parallel plate system for various gaps indicated that the plate's roughness did suffice 
to prevent the wall slip (data not shown here). 
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Figure 2: (a) Schematic view of the flow control system: TS - vertical translational stage, M - stepping motor, GB - gear box, GB - micro step 
controller, PC - personal computer, IFC - inlet fluid container, OFC - outlet fluid container, FC - flow channel, FC - flow channel, Ii - working 
fluid inlet, Oi - working fluid outlet, B - digital balance, (b) Schematical representation of the hydrostatic pressure ramp, is the characteristic 
forcing time (the averaging time per stress point) and N is the total number of steps. 



A second concern was related to the possible artefacts introduced by fluid evaporation during long experimental 
runs. In order to prevent this, a solvent trap has been placed around the free fluid meniscus. The sealing of the solvent 
trap on the base plate of the rheometer has been insured by a thin layer of vacuum grease. After each experimental 
run it has been carefully checked (by visual inspection) that no significant changes in the shape of the meniscus 
occurred. Additionally, we have checked at the end of each run that one can reproduce the viscosity measured during 
the pre-shear step which indicated us that the evaporation effects were either minimal or absent. 

2.4. Flow control 

A slow controlled pressure flow has been generated by controlling the hydrostatic pressure difference between the 
inlet and the outlet of the flow channel as schematically illustrated in Fig. p_| The working fluid is held in a container 
IFC with a large free surface rigidly mounted on a vertical translational stage TS. The translational stage is driven by 
a step motor M via a gear box GB. 

The position and the speed of the stepping motor are controlled by a micro-step controller (model IT 116 Flash 
from Isel Gmbh) connected to a personal computer PC via a RS 232 serial interface. To control the motion of the 
stepping motor, a graphical user interface (GUI) has been developed under Matlab® . The vertical position of the inlet 
fluid container IFC is controlled with an accuracy of 37 jim. After passing through the channel, the fluid is collected 
in a outlet fluid container OFC placed on a digital balance B. The balance provides real time reading of the mass of 
the fluid discharged through the channel via a RS 232 serial port interfaced under Matlab®. 

The pressure driving the pipe flow is p(f) = pgh{t) — p t where p = 1050 kg m -3 is the density of the Carbopol® 
solution, g the acceleration of gravity and p, the pressure drop in the connecting tubes. Measurements of the flow rate 
obtained from the mass of the discharged fluid (the data are not shown here) allowed one to estimate that the pressure 
drop along the connecting tubes p, roughly accounts for 5% of the total hydrostatic pressure difference. 

The imposed hydrostatic pressure mimics the flow ramps employed in Ref. ^utz and Burghelea (2009) to investi- 



gate the yielding phenomenon and consisted of linear hydrostatic pressure ramps for both increasing and decreasing 
values of the height difference h(t) (see Fig. ^|(b)): 

h(t) = £ C[i] (H(t - it ) - H(t - (/ + l)/ )) (1) 

;=o 

where 

C\i\ = , ma * I 1 ~ 2 (2) 

u (N-l)t [N-{i+l), o.w. 

and H(t ) is the Heaviside function. 

The pressure ramp described by Eq. p] and illustrated in Fig. p| (b) consists of N steps each lasting a time fo. As 



in the case of the stress ramps studied in Putz and Burghelea (2009) fo will be referred to as the characteristic forcing 
time. 
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2.5. Flow field measurements 

Time series of the velocity fields were obtained by a iterative multi-grid Digital Particle Image Velocimetry (DPIV) 
algorithm implemented in the house under Matlab® Scaranoand Rhiethmuller (2001); Raffel etal. (September 2007). 
For this purpose, simultaneously with the hydrostatic pressure ramp illustrated in Fig. |2j(b) a sequence of flow images 
has been acquired at a frequency of 10 frames per second (fps). Corresponding to low values of the driving pressure, 
several flow images have been skipped in order to increase the inter-frame and maintain the average displacement of 
the tracer particles in the optimal range of 5 to 15 pixels, Raffel et al. (September 2007); Scarano and Rhiethmullei 
(2001 ). Together with this, the smallest interrogation window has been adapted to the mean flow velocity correspond- 
ing to each step of the pressure ramp. The spatial resolution of the measured flow fields was 99.5 Um. To account 
for the significant mean flow speed variation during the controlled pressure ramps illustrated in Fig. 0(b), an adaptive 
inter-frame between the subsequent images passed to the DPIV algorithm has been used. Thus, as the pressure drop 
has been increased, the inter-frame was gradually decreased from 1.5 s (corresponding to the lowest value of the 
pressure drop) to 0.1 s (corresponding to the largest value of the pressure drop investigated). This procedure allowed 
one to maintain the instrumental errors in the desired confidence interval, regardless the mean flow speed. 

The accuracy of the method has been carefully checked by running test measurements for low Reynolds number 
Poiseuille flows with Newtonian fluids and comparing them with the analytical solution. Thus, it has been checked 
that using this adaptive DPIV protocol the instrumental error of the measured velocity fields does not exceed 7% 
within the entire range of driving pressures we have explored. 



3. Experimental results 

3.1. The solid-fluid transition in a rheometric flow: the role of wall slip 

Prior to characterising the unsteady pipe flow of a Carbopol® gel we focus on the solid-fluid transition which we 
have investigated in detail in |Putz and Burghelea (|2009| ) for this material. It is known that Carbopol® gels exhibit 



slip while flowing past smooth surfaces, |Piau ( 2007| ). Thus, as the flow channel used in this study is glass-made it 



is highly expected to observe a non-negligible velocity slip at the wall. This motivated us to perform a comparative 
rheological investigation of the Carbopol® gel in two separate geometries with and without wall slip. The rheological 
response of the Carbopol® sample is tested during a stepped stress ramp performed at a constant temperature. This 
stress ramps closely mimics the pressure ramp that drives the channel flow illustrated in Fig. ^(b). The rationale of the 
choice of this rheological protocol is to be able to correlate the rheological measurements with the flow measurements 
performed in the channel which follow a similar forcing scheme. 




Figure 3: (a» Flow curves measured with a polished geometry (□, ■) and with a rough geometry (o, The full/empty symbols refer to the 
increasingmecreasing branch of the flow curves. The full line is a Herschel-Bulkley fit that gives T v = 0.64 ± 0.003(Pci), K = 0.4 ± 0.002(P«j"), 
n = 0.54 ± 0.001. ibl Dependence of the hysteresis area on the characteristic forcing time ?o measured with a smooth geometry (□) and a rough 
one (■). The full line (— ) is a guide for the eye, P °c tZ 63 and the dash dotted line (— . — ) is a guide for the eye, P °= ^°'° 3 - 
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The results of these measurements are presented in Fig. |[ In Fig 
measured with both rough plates (the rhumbs, (o, ♦)) and smooth plates (t 



3(a) are displayed controlled stress ramps 
))■ 



the squares, (□, I 

As already emphasised in |Putz and Burghelea| (2009); Weber et al. (2012), our choice of representing the Theolog- 
ical data in the coordinates y = /(t) (which might seem somewhat unusual to some) is motivated by the fact that this 
representation allows a more accurate identification of the various deformation regimes and a rigorous assessment of 
the reversibility of the deformation states. 

For the case of the measurements in the absence of all slip (the rhumbs (o, ♦) presented in Fig. 3(a)) and cor- 
responding to low values of the applied stress ( T < 0.64Pa ) the measured rate of deformation is constant for both 
branches (increasing/decreasing) of the stress ramp and, consequently, the strain 7 is a linear function of the applied 
stresses which is the signature of a elastic solid deformation regime []. Within this regime the deformation states are 
not reversible upon increasing/decreasing applied stresses which translates into different elastic moduli: G u > Gd- 
Here the indices denote the increasing/decreasing stress ramps, respectively. 

For large values of the applied stresses ( T > 0.64-Pa) a reversible fluid regime is observed. One has to emphasise 
that the transition from a irreversible solid like deformation regime to a reversible fluid one is not direct, but mediated 
by an intermediate deformation regime which can not be associated with neither a solid like behaviour nor a fluid one 
but with a coexistence of the two p hases. A similar smoo t h soli d-fluid transition has been recently observed in a low 
Reynolds number steady pipe flow, ^erez-Gonzalez et"aL ( 2012 ). 

The irreversibility of the deformati on states observed wi t hin th e solid and the intermediate deformation regimes 
and i nitially repor t ed by Putz et al in Putz and Burghelea| ( 200S ) has been recently confirmed by others, Divoux 
et al. ( ]2010l [20 1 1[ ); |Divoux et al.| ( |2012| ). A similar irreversible rheological behaviour has been found in rheological 
studies perfor med on various grades o f Carb o pol® gels us i ng var ious rheometers (and various geometries) at various 
temperatures, Moyers-Gonzalez et al. (2011); 



Weber et al 



(2012). 



We note that, during ea c h of the rheological measureme nt s reported in t his m anuscript and in our previous studies 
^utz and Burghele"a ( 2009 ); Moyers-Gonzalez et al. ( 2011 ); Weber et al. ( 2012 ), the Reynolds number was s maller 
than unity, therefore the experimentally observed irreversibility of the deformation states illustrated in Fig. 3(a) should 
not be associated with inertial effects. 

Thus, although Carbopol® gels have been long considered as "model " yield stress fluids that display no thixotropic 
effects and their rheological behaviour can be accurately described by the Herschel-Bulkley picture, this recently 
found irreversible rheological behaviour should be accepted as a physical fact and consequently trigger new devel- 
opments that can properly account for this picture and its impact on the description of various realistic viscoplastic 
flows. 

To understand the role of wall slip in the yielding of the Carbopol®, we have performed the same type of rheolog- 
ical measurements with smooth parallel plates. The main effect of the wall slip is to shift the solid-fluid coexistence 
regime to lower values of the applied stresses (see the squares (□, ■) in Fig. 3(a)). The apparent yield stress measured 
in the presence of wall slip t" is nearly an order of magnitude smaller than the yield stress measured in the absence 
of slip, Zy. The data acquired with and with out wall slip overlap only within the fluid regime, corresponding to shear 
rates 7 > y c (with % rj 0.8i _1 see Fig. 3(a)). This critical value of the shear rate above which the slip and no slip data 
overlap is comparable to th e critical shear rate f ound by Bertola and his coworkers below which no steady state flow 
could be observed by MRI, Bertola et al. ( 2003 ). The existence of the critical shear rate % is also consistent with the 
measurements of the velocity pro files performed by S almon and his coworkers with a concentrated emulsion sheared 
in a small gap Couette geometry, Salmon et al. ( 2003 ). 

Additionally, we note that it is solely within the viscous d eform ation regime that all the data sets can be reli- 
ably fitted by the Herschel-Bulkley model (the full line in Fig. f3(a)|) which reinforces the idea that only within this 
deformation range the Carbopol® gel behaves like a simple or "model" yield stress fluid. 

A next fundamental question that deserves being addressed is how the irreversible character of the rheological 
flow curves observed within the solid-fluid coexistence regime in the form of a hysteresis loop is related to the rate 
at which the deformation energy is transferred to the material or, in other words, to the characteristic forcing time 
fo. It is equally important to understand if the presence of wall slip influences this dependence. To address this 



'Note that the stress ramp is linear in time as illustrated in Fig. g|and, therefore, one can conclude that the material deforms as an elastic solid 
within this regime, G = T/y = ct 
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points we compare measurements of the area of the hysteresis observed in the flow curves presented in Fig. 3(a) 
P = J y"dAz" — J y^if At'' performed for various values of to for both slip and no-slip cases. Here the indices "u,d 
refer to the increasing/decreasing stress branches of the flow ramp. As already discussed in Ref. f*utz and Burghelea 
( 2QQ9| ) the area P has the dimensions of a deformation pow er de ficit per unit volume of sheared material. The results 

In the absence of wall slip , the deformation power 
(the full squares (■) in Fig. 3(b)). A similar scaling 



P(b) 



of these comparative measurements are presented in Fig. 
deficit scales with the characteristic forcing time asP=< t 
has been found in the absence of wall slip for several Carbopol® solutions with various concentrations and at various 
temperatures in Ref. |Putz and Burghelea ( 2009 ),F«f n '. Th e difference in the scaling exponent found in the present 



Putz and Burghelea] ( |2009 ) may be related to differences in the grade of the 



study and the one initially reported Ref. 
Carbopol®. 

The existence of a rheological hysteresis for Carbopol® gels characterised by a decrease of the power deficit with 
the characteristic forcing time cons istent with both our initial finding reported in Ref. Putz and Burghelea (2009) 



and the data presented in Fig. 3(b) was afterwards confirmed by others i n independent experi ments performed with 



Divoux et al 



(2012) 



different Carbopol® solutions and using different experimental protocols, 

The next question we address is what is the influence of the wall slip phenomenon on the rheological hysteresis 
observed during stepped stress ramps and on its scaling with the characteristic forcing time to- To address this points 
similar measurements of the hyste resis area P for various values of to have been performed in the presence of wall 
slip (the empty squares (□) in Fig. 3(b)). The presence of the wall slip affects both the magnitude of the deformation 



power deficit and its scaling with the characteristic forcing time by nearly suppressing it: P t Q (the dash-dotted 



line(- -)inFig. 3(b)) 



The central conclusion of this comparative investigation of the solid-fluid transition in a rheometric flow in the 
presence and in the absence of wall slip is two-fold. Firstly, the wall slip phenomenon and the irreversible deformation 
are observed in the same range of applied stresses and measured rates of shear. This means that, for this particular 
flow, the irreversible solid-fluid transition and the wall slip are coupled phenomena. Secondly and reinforcing the first 
part of the conclusion, the emergence of the wall slip modifies the irreversible solid-fluid transition by suppressing the 
dependence of the hysteresis power loss on the rate at which the material is being forced, 1 /to- 

From a theoretical point of view, the experimentally observed coupling between the wall slip and the yielding 
phenomena in a unsteady rheometric flow with slippery boundaries is still insufficiently documented and complicates 
the physical picture. 

I n the absence of wall slip, a simple model able to describe the irreversibility of deformation states observed in 
has been recently propos ed, ^utz and Burghelea (2009). As previously suggested by several authors (Moller 



Fig. 
et al 



3(a) 



I2TJ06); 



Dullaert and Mewis 



(2006)), the yilding process of a physical gel sample under shear was interpreted in 
terms of a dissociation reaction, S ^ S + F and be modelled by the following kinetic equation: 



d<i> 

— =R d (®,t,x)+R r (®,t,'c) 



(3) 



where S, F denote the solid and fluid phases, respectively, 4> = [S] is the concentration of the solid phase, T is the 



forcing parameter (the applied stress), Rj = —K c i 



1 + tanh 



( T ~ T y 
\ w 



<t> is the rate of destruction of solid units and 



1 — tanh f^jnj (1 ^ < t > ) is the rate of fluid recombination of fluid elements into a gelled structure, Putz and 



R r — K, 
Burghelea (|2009j). 

As a constitutive equation a thixo-elastic Maxwell (TEM) type model was used: 



G dt 



T= 77(7)7 



where the viscosity of the fluid phase is given by a regularized Herschel-Bulkley model, rj (7) = K | y\" 



-1 



(4) 

1— exp(— m\y\) 



y 1 71 

Here G is the static elastic modulus, K the consistency, n the power law index and m is the regularisation parameter. 

The choice of this constitutive equation is motivated by the pres ence of elastic effects in the intermediate defor- 
mation regime (see the cusp in decreasing stress branch in Fig. 3(a) and the corresponding discussion). We note that 
in the limit <t> — > 1 equation ^reduces to Hooke's law, T = Gy (which describes the solid deformation regime), and in 
the limit <t> — > it reduces to a regularized Herschel-Bulkley model (which describes the viscous regime). 
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In spite of its simplicity, the model briefly presented above proved its ability of describing a large number of 
rheological flow curves similar to that presented in Fig. 3(a) (the squares (□, ■)) measured in the absence of wall slip 
for Carbopol® solutions of various grades, concentrations and at various temperatures, ^utz and Burghele"a ( |20Q9| ): 
Weber et aT] ( [2012 ). Additionally and quite remarkably, for a given aqueous gel solution (0.1%, Carbopol® 940) and 
for the same choice of the model parameters (t v , G, n, K u , Kj, w), the model accurately fits both controlled stress 
increasing/decreasing ramps and controlled stress oscillatory sweeps in a large range of oscillation amplitudes and 
frequencies, Putz and Burghelea ( 2009| ). 

The presence of wall slip whi ch art ificially diminishes the apparent yield stress t" < T v and significantly alters the 
shape of the flow curves (see Fig. 3(a)) renders the simple model presented above inapplicable in its present form as it 
contains no terms that can account tor the slip behaviour and for its coupling to the solid-fluid transition. A modified 
version of this model that can account for wall slip by introducing a delay term in the phase equation Eq. [3] will be 
described elsewhere. 



3.2. The solid-fluid transition in a unsteady pipe flow 

Following the comparative investigation of the yielding of a Carbopol® gel in a rheometric flow in the presence 
and in the absence of wall slip presented in Sec. pTl| , the main question that arises is to what extend the findings 
on the solid-fluid transition investigated in this particular laboratory scale flow could be transferred to flows that are 
more relevant from a practical perspective, such as a slow (inertia-free) unsteady pipe flow. Of equal importance is 
understanding the role of wall slip in the dynamics of the solid-fluid transition. 



3.2.1. Dynamics and structure of the unsteady flow fields 

Within the framework of the Herschel-Bulkley model and in the absence of wall slip, the mean flow velocity U a 
is related to the driving pressure drop via: 




! + !/« 



1 



S(i+V») 



S<1. 
S>1. 



see e.g. Bittleston et al. ( 2002 ). The dimensionless parameter t, is defined: 

^ ILXy 



(5) 



(6) 



The Eq. g relies heavily on the Herschel-Bulkley yielding picture, the flow steadiness and the absence of velocity 
slip at the channel's wall. Although we expect that none of these assumptions are fully fulfilled by the flow we 
investigate, we still believe it is reasonable to use it in order to get just an indicator for order of magnitude for the 
applied pressure drop that corresponds to yielded states, Ap y , which would ultimately indicate us a suitable range of 
the driving pressure drops Ap that needs to be explored during our experiments in order to systematically characterise 
the solid-fluid transition for this type of unsteady pipe flow. 



Using the yielding criterion ^ w 1 for case of the flow channel described in Sec. 2.1 and the Carbopol® solution 
characterised in Sees. 2.2, 2.3, the pressure drop corresponding to the yielded (flowing) states can be estimated as 



Ap y w 774Pa. 

To characterise the solid-fluid transition in the flow channel, a time series of flow fields has been acquired during 
a controlled increasing/decreasing pressure ramp (see Fig. §(b)). This controlled pressure ramp is similar to the 
controlled stress ramps used to characterise the solid-fluid transition in a rheometric flow. 

The corresponding time series of the absolute value of the centre line velocity \U p \ is presented in Fig. The 
time interval during each the applied pressure was maintained constant (the characteristic forcing time fo) was fixed 
at 7.5 s. 

In order to capture the details of the solid-fluid transition in an unsteadily forced flow, the maximal value of the 
applied pressure drop has been chosen nearly three times larger than the above estimated pressure drop corresponding 
to the yielded state, Ap y . 
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Figure 4: Time series of the absolute value of the plug velocity U P = U . The insert presents a magnified view of the time series after the recoil 
effect (see the text for the explanation) is observed. The characteristic forcing time was to = 1.5s. 



The spikes visible in the velocity time series presented in Fig. ^ correspond to the change of the driving pressure 
realised by the displacement of the vertical translational stage (TS in Fig. ||(a)) which carries the inlet fluid container 
IFC. Its magnitude depends on the inertia of the mechanical system and the acceleration of the stepping motor M. 

To avoid accounting for this response of the mechanical system which controls the pressure drop along the flow 
channel, the flow data acquired during these transients (extended over roughly 0.5s in each of the experiments reported 
here) has been discarded. 

Corresponding to the first five values of the driving pressure on the increasing branch of the ramp the time depen- 
dence of the centreline velocity does not reach a steady state during the characteristic forcing time to = 7.5s. 

A further increase of the driving pressure reveals an interesting feature of the velocity time series in the form of 
a non-monotone time dependence: upon the change of the pressure drop, the centreline velocity first decreases and 
then increases. A phenomenological explanation of this observation may be given in terms of a strong heterogeneity 
of the flow, i.e. a dynamic coexistence and competition between yielding of solid material elements and the re- 
combination of yielded fluid elements which, as illustrated by the rheometric data presented in Fig. 3(a), occurs 
within an intermediate range of the applied stresses (pressure drops). Indeed, as the driving pressure drops Ap are 
further increased, this non-monotone behaviour is no longer observed and the plug velocity reaches a steady state 
during a characteristic time smaller than to- This indicates that, within this range of driving pressure, a steady yielded 
state is achieved, Fig. |j. 

The time series presented in Fig. || may also provide a first indication on the reversibility of the flow states upon 
increasing/decreasing the driving pressure drops. Thus, one can note that the data is not symmetric with respect to 
the vertical line t /to = 10. This observation is consistent with the irreversibility of the deformation states observed 
for the rheological data acquired at low and intermediate values of the applied stresses, Fig. 3(a). An even clearer 



signature of the irreversibility of the deformation states can be noted if one monitors the time series data acquired 
during the last four steps of the decreasing branch of the pressure ramp and highlighted in the insert of Fig. [|. For 
these steps of the pressure ramp velocity fluctuates strongly and, corresponding to the the last step a flow reversal is 
apparent: the absolute value of velocity seems to increase as the pressure drop is decreased to its smallest value. A 
careful inspection of the flow images acquired during this step clearly indicate a reversal of the flow direction. This 
rather unexpected effect can be associated to a elastic recoil effect similar to the recoil obs erved during the rheological 



measurements on the decreasing stress branch and manifested by the cusp visible in Fig. 3(a) 



We focus in the following on the structure of the unsteady flow fields and their fluctuations. Profiles of the absolute 
value of the mean flow velocity U av and the reduced velocity fluctuation U rms /U av measured for several values of the 
driving pressure Ajj on both the increasing and the decreasing branch of the flow ramp illustrated in Fig. ^|(b) are 
presented in Fig. 
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r/R r/R 

Figure 5: (a) Profiles of the absolute value of the mean flow velocity U av . (b) Profiles of the reduced root mean square variation of the velocity, 
Unm/Uay. The symbols refer to different values of the pressure drop: (□,■)- Ap = 310.8Pa, (o, .) - Ap = 621.6Pa, (A, A) - Ap = 932.3Pa . The 
full/empty symbols refer to the data acquired on the increasing/decreasing branch of the pressure ramp, respectively. The characteristic forcing 
time is to = 10s. 



For each value of the driving pressure, on the increasing branch of the flow ramp the mean velocity profile indicates 
the presence of a central unyielded plug, Fig. [| as typically expected for a viscoplastic fluid. Also, for each value of 
the pressure drop, a slip behaviour is observed near the channel walls, r = ±1. For the largest values of the driving 
pressure (Ap = 932. 3Pa and Ap = 621. 6Pa, (the triangles (A, A) in Fig. || (a))) the velocity profiles measured on 
the increasing/decreasing branches of the ramp coincide consistently with the time series measurements presented in 
Fig. Q This indicates that a stationary fluid flow regime has been achieved. Within this state the velocity fluctuations 
quantified by U rms /U av are the smallest and mainly localised near the channel wall which is another indicator for the 
steadiness of the flow (see Fig. || (b)). 

Corresponding to the lowest driving pressure, however, the profile measured on the decreasing branch (the empty 
squares, □, in Fig. ^ (a)) falls clearly below the profile measured on the increasing pressure branch. This indicates 
that within this range of pressure drops the flow is irreversible upon increasing/decreasing the forcing. It is equally 
worth noting that, corresponding to this driving pressure, the mean flow profile is less smooth around the centreline 
of the channel which indicates a rather large level of velocity fluctuations. It also indicates that smaller the driving 
pressure is further the flow is from a steady state. 

To probe this hypothesis we turn our attention to the transversal profiles of the reduced velocity fluctuations 
Urms/Uav, Fig. ^| (b)). For large driving pressure (the triangles in Fig. || (b)) the reduced velocity variations account 
for roughly 6% of the mean flow velocity in the bulk region of the flow (around the unyielded plug) and they may reach 
up to 20% of the mean flow velocity in the wall region (due to the smallness of the measured velocity). The rather 
large values of the fluctuations observed near the channel walls are instrumental and they are due to the smallness of 
the measured velocity within this regions (here, the DPIV heavily relies on the sub-pixel interpolation because of the 
very small displacements of the flow tracers). 

On the decreasing branch of the controlled pressure ramp and for low driving pressures a localisation of the 
velocity fluctuations can be observed in the form of a local maximum centred around the middle of the channel 
(r/R = 0) can be observed (the empty circles (o), the empty triangles (A) and the empty squares (□) in Fig. §(b)). 
The localisation of the velocity fluctuations near the centre line of the flow channel may be understood in terms of 
large fluctuations of a elastic (un-yielded) plug which is consistent with the experimental findings for the case of a 



oscillatory pipe flow, Park and Liu (2010) 



3.2.2. Reversibility of the flow states, the role of wall slip 

We focus in the following on the reversibility of the flow states upon increasing/decreasing pressure drop and its 
connection to the irreversible solid-fluid transition observed in a rheometric flow in the presence of wall slip illustrated 
in Fig. § 
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Figure 6: Dependence of the absolute value of the mean flow velocity U av on the applied pressure drop Ap for various values of the characteristic 
forcing time: (a) to = 4s, (b) to = 5s, (c) to = 1.5s, (d) = 10s, (e) to = 15*. 0>) = 20s. The squares (□) refer to the increasing branch of the 
stress ramp and the circles (o) refer to the decreasing branch of the stress ramp. 



We note that our flow investigation technique required us to use a flow channel with slippery walls so a comparison 
of the slip and no slip cases cannot be performed for this pressure driven flow and, consequently, the yielding transition 
can not be decoupled from the wall slip. 

Such a comparative study remains open to future experimental investigations that use measuring techniques that 
do not require optical transparency of the flow channels such as ultrasonic velocimetry. 

Measurements of the absolute value of the mean flow velocity U av performed for increasing/decreasing values of 
the driving pressure Ap and several values of the characteristic forcing time are presented in Fig. ^. The mean flow ve- 
locity U av has been calculated by numerical integration of the DPIV measured velocity profiles, U av = 2n L rU[f)dr 

B 

Regardless the value of the characteristic forcing time, the data acquired for increasing/decreasing driving pres- 
sures overlap only within the fluid regime Ap > Ap y . This indicates that only corresponding to this range of the driving 
pressures the Carbopol® gel is fully yielded. Below this critical value of the driving pressure a hysteresis of the flow 
states is observed. This observation is fully consistent with the rheological hysteresis observed in the presence of wall 



slip (the squares (□) in Fig. 3(a)). As previously mentioned through our paper, the cusp visible on the decreasing 



pressure branch corresponds to a elastic recoil effect manifested by a reversal of the flow direction. This indicates 



The axial symmetry of the flow has been assumed. 
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Figure 7: Dependence of the area of the hysteresis of the dependencies of the mean flow velocity on the driving pressure drop presented in Fig. 
on the characteristic forcing time fo. The full line is a guide for the eye, A <x tg° 1 . 



that, as in the case of the solid-fluid transition observed in a rheometric flow f'utz and Burghelea (2009), the elasticity 
cannot be ignored while studying the unsteady yielding in a slow pipe flow flow of a Carbopol® gel. 

The dependence of the hysteresis area of the dependence of the mean flow velocity U av on the driving pressure 
Ap defined by A = jU^dAp" -JU% v dAp d is illustrated in Fig. | is presented in Fig. \A Here the indices "u, d" 
refer to the increasing/ decreasing branch of the controlled pres sure r amp sketched in Fig. g (b). As in the case of the 
rheological hysteresis in the presence of slip illustrated in Fig. |3(b)| a weak dependence of the hysteresis area on the 
characteristic forcing time is observed, A « fcT 0-1 . 

This fact indicates that the presence of wall slip modifies the irreversibility of the flow curves within the solid-fluid 
coexistence regime and prompts us to focus closer on the slip behaviour observed during the stepped pressure ramp 
and its relation to the rheological behaviour illustrated in Fig. B. 

The dependence of the slip velocity U s on the driving pressure drop Ap measured for several values of the charac- 
teristic forcing time fo is presented in Fig. || 

As for the pressure dependence of the mean flow velocity U av illustrated in Fig. || and discussed above, one can 
clearly see in Fig. || that, regardless the value of fo, the dependence of the slip velocity on the applied pressure drop 
Ap is not reversible upon increasing/decreasing pressure over the entire range of pressures. This is another indicator 
that the slip behaviour can not be entirely decoupled from the irreversible solid-fluid transition observed during the 
rheological measurements presented in Fig. 3(a). Regardless the value of the characteristic forcing time fo, within the 
fluid flow regime (Ap > 500Pa) the slip velocity scales with the applied pressure drop as U s °< Ap 2 (see the full lines 
in Fig. ||). This scaling law differs from the one found by Gonzalez and his coworkers in Perez-Gonzalez et al. ( 20 1 2 ). 
U s « Ap 876 (see Fig. 6 in Ref. Perez-Gonzalez et al. (2012) in connection to their Eq. (3)). As the scaling law we 
found is practically insensitive to the characteristic forcing time, we believe that this discrepancy does not originate in 
the steadiness of the pipe flow investigated in Perez-Gonzalez et al. (2012) but in the different rheological properties 
of the Carbopol® solutions they have used (nearly an order of magnitude difference in both the yield stress and the 
consistency). Thus, from this comparison, it is apparent that even within the reversible fluid regime (t > T v ) the wall 
slip behaviour remains correlated to the rheological properties of the solution. 

To get a deeper insight into this correlation, we focus on the dependence of the slip velocity U s on the wall 



velocity gradients ^ 



, Fig. H The measurements presented in Fig. || are performed for various values of the 

>=±« n 
characteristic forcing time fo on both the increasing and the decreasing branch of the pressure ramp (see Fig. |2J (b)). 

For each of the data sets presented, the wall velocity gradients have been obtained by numerical differentiation of 

the DPIV measured averaged velocity profiles near the wall and the slip velocity has been obtained by extrapolating 

the same time averaged profiles at the wall (see Fig. || (a)). To increase the accuracy of the numerical differentiation, 

a central difference differentiation scheme has been used. 
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Figure 8: Dependence of the slip velocity U s on the applied pressure drop Ap for various values of the characteristic forcing time: (a) to = 4s, (b) 
to = 5s, (c) to = 7.5s, (d) fo = 10s, (e) to = 15s, (b) fo = 20s. The squares (□) refer to the increasing branch of the stress ramp and the circles (o) 
refer to the decreasing branch of the stress ramp. The full lines in each panel are guides for the eye, U s « Ap 2 . 
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Figure 10: Dependence of the slip velocity U s on the wall shear stresses T w . The full line is a power law fit, U s = 2.3 ■ 10 (±5.4- 10 



The dotted and the dash-dotted lines are the theoretical predictions by Piau (see Ref. 3 iau (2007)): U s « T„ 
t = 4 (s), (o, .) - to = 7.5 (j), (□, ■) - to = 10 (j), (V, A) - / = 15 (s), (<,-*)- (o = W(sJ. 



1/3 



t/ s °= T^. The symbols are: (A, A) ■ 



Previous w orks have reported a decrea se of the slip velocity with the applied stresses (velocity gradients) beyond 
the yield-point, Nguyen and Bogei ( 1992 ). Contrarily to this, our experimental findings indicate a monotone increase 
of the slip velocity with the wall velocity gradients, regardless the deformation regime (solid, fluid or solid-fluid) and 
the degree of flow steadiness (the value of fp). This con c lusion is consistent, however, with the experimental findings 
for the case of a steady pipe flow, Perez-Gonzalez et al. ( 2012 ). 

Regardless the deformation regime, the value of the characteristic forcing time to and the type of the pressure 
ramp (increasing or decreasing pressures) a universal dependence of the slip velocity U s on the wall velocity gradients 

is observed, Fig. ^, and all the data follow a line who's slope indicates a slip length of roughly 300^im (the 

r=±R 

full line in Fig. 0) 

Thus, one can conclude from theses measurements that the dependence of the slip velocity on the wall velocity 
gradients has an universal character and uniquely defines a slip length which is independent on the deformation 
regime but may depend on the roughness of the channel wall and on the physicochemical properties of the Carbopol® 
solution. Two natural but fundamental questions arise at this point: 

1. What is the extent of this universal behaviour? 

2. Is there a simple phenomenological picture underlying the wall slip behaviour observed for a Carbopol® gel? 

A first phenomenological picture for the wall slip phenomenon (observed in suspensions) has been suggested in 
early 20' s by Bingham by the following quote: "Slip comes from a lack of adhesion between the material and the 
shearing surface. The result is that there is a layer of liquid between the shearing surface and the main body of the 
suspension" , see page 231 in Ref. Bingham (1922). 

Kaylon developed further on the slip picture suggested by Bingham and proposed the following relationship 
between the slip velocity and the wall shear stress, Kalyorj ( [2005 ): 



K=j3(TW) 1/ns (7) 

The slip coefficient j5 is related to the width 8 of the slip layer, the consistency K s and power law index n s of the 
interstitial fluid arrested within the slip layer (referred to in Ref. Kalyon| (2005) as the "binder fluid") via /3 = j^. 

Carbopol® micro gels are not suspensions in the strict sense of the term but they can still be viewed as a jammed 
system of swollen gel micro-particles, Piau (2007). Thus, bearing in mind that the gel was prepared in an aqueous 
solvent, it looks reasonable to us to consider the "binder fluid" as Newtonian, K s = 1 mPas,n s = 1. 

To test the applicability of Eq. ^ with the assumptions above for our data, we re-plot the slip velocity data presented 



in Fig. ^ versus the wall shear stress defined by T w 



dv 
777 



r=±R 
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r=±R 



For this purpose we have related the 



wall velocity gradients measured by numerical differentiation of the DPIV measured velocity profiles to the wall shear 
stresses via the rheological flow curve measured in the presence of wall slip and presented in Fig. 3(a) (the squares) 
with the wall velocity gradient playing the role of the rate of shear y. The result of this evaluation is presented in Fig. 
|To[ Within this representation the universality observed in Fig. ||is broken and three different slip regimes can now be 
clearly observed. Corresponding to values of the wall shear stress signi ficantly smaller then the apparent yield stress 
measured in the presence of wall slip (x w < t" w Q.lPa, see Fig. 3(a)) the measured slip velocity displays a rather 
large scatter, most probably related to the smallness of both the slip velocity and wall velocity gradients. However, in 
spite of this data scatter and as wall shear stress approaches the apparent yield stress t" measured in the presence of 
wall slip, one can note a non-monotone dependence of the slip velocity on the wall shear stress w hich is a reminiscent 
of the cusp observed in the dependence of the rate of shear on the applied stress presented in Fig. 3(a). This behaviour 
clear ly departs from a Navier-like slip behaviour and it is related to the non-trivial yielding process illustrated in Fig. 
3(a) which occurs via intermediate range of stresses where the solid and fluid phases coexist, |Putz and Burghelea 
( 2009 ). Above the apparent yield stress and corresponding to the fluid deformation regime a nearly linear scaling 
of the slip velocity with the wall shear stress is observed, U s « T ^-32±0.26 ( me f u jj jj ne m pjg indicates 
a Navier-like slip behaviour which can be described by Eq. ^. This nearly linear scaling reinforces the assumption 
above that the "binder fluid" is Netwonian, i.e. n s « 1. We note that a similar scal i ng law has been found in the recent 
micro-flow experiments by Gonzalez and his coworkers, |Perez-Gonzalez et al. ( 2012 ) (see Fig. 6 in their paper): 



oc T 



0.876 



The pre-factor of this linear dependence is j3 ~ 0.00023 mPa s which leads to an estimate for the 
thickness of the Newtonian shear layer 8 « 0.23/im. This estimate is fairly close to the value measured by Jiang and 
his coworkers for the case of hydroxypropyl guar (HPG) gel, 8 « 0. \\lm. 

It is interesting to compare our experim ental l y foun d scaling of the slip velocity with the theoretical prediction of 
Piau for a Carbopol® gel presented in Ref Piau ( 2007 ). Using a Hertzian framework to model the contacts between 
neighbouring m icro-gel elastic particles (a detailed description of the theoretical framework can be found in Ref. 
Johnson (1985)), Piau predicts two scaling laws of the slip velocity with the wall shear stress depending on the 

9 I /3 

packing of the soft particles: U s « T„, for a closely packed system and U s V in the case of a loosely packed 
system. The scaling exponent obtained by fitting the data presented in Fig. [Io| corresponding to the fully yielded 
regime differs from these theoretical predictions (see the dotted and the dash-dotted lines in Fig. [Io| ). 

A possible reason of this discrepancy might originate in the fact that the interstitial lubricating fluid considered 
m Ref. piau| ( p007| ) was non-Newtonian and, to get a deeper insight into the slip behaviour of Carbopol® gels a 
additional discussion of this issue is in order. A physical rationale for the choice of a non-Newtonian (shear thinning) 
interstitial fluid in Ref. Piau (2007) might be to decouple at a micro scale the shear thinning behaviour (by attributing 
it to the interstitial fluid) from the yield stress behaviour by implicitly attributing it to the elastic-solid Carbopol® 
micro particles. 

The disagreement we pointed out above seems to indicate that this decoupling is not as straightforward and should 
be more cautiously considered in future attempts to theoretically describe the relation between flow properties, wall 
slip and microstructure. Future theoretical approaches should definitely reconsider this problem. 

A more elaborate d theory of the slip dynamics observed in a soft gel was presented by Meeker and his coworkers, 
[Meeker et al. ( 2004 ). Due to the non-negligible level of scatter of the data presented in Fig. induced by the 
numerical differentiation of the measured velocity profiles, however, a direct comparison of our results with these 
predictions is not possible. 

Currently, there exist no experimental studies of the microscopic scale dynamics of Carbopol® gels under shear 
in the vicinity of slippery surfaces able to confirm the existence of a shear lubricating layer and provide a quantitative 
measure of its thickness. Thus, the estimate of the width 8 of the lubrication layer we have obtained indirectly from 
the analysis slip velocity data should only be regarded as plausible. A direct experimental investigation by fluorescent 
visualization of the gel network as recently performed by Gutowski et al Gutowski et al. (2012) would be highly 
needed to test this estimates. We can, however, provide some additional arguments for the meaningfulness of this 
estimate as follows. 

Systematic microfluidic investigations of the diffusion of nano-tracers within a Carbopol® micro-gel indicate that 
the average characteristic size of gel micro particles is of the order of R ar w Iflm, pppong et al. (2006); Oppong and 
de Bruyn (2007). A width of the shear layer 8 f=s Q.23R av together with the assumption of a random close packing of 
the micro-gel particles is consistent with a rough standard deviation of the particle size R st d/Rav ~ 23%. Again, based 
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on the microscopic experimental evidence presented in Refs. Oppong et al. (2006); Oppong and de Bruyn (2007) such 
a value can be considered, at least for now, plausible. 

To conclude this part, our experimental data seems to confirm that in a fluid state (t > t") the Bingham inspired 
phenomenological picture of a slip layer proposed by Kaylon is applicable to the unsteady pipe flow of the Carbopol® 
gel investigated through our paper. 

Another important fact that needs to be pointed out regarding the data presented in Fig. 
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is that, regardless the 

value of the characteristic forcing time to and for both increasing/decreasing values of the driving pressure Ap (the 



full/empty symbols in Fig. ttOl respectively), the slip velocity data overlap onto a single master curve. This invariance 
of the slip velocity data with the characteristic forcing time to can be explained b y the weak dependence on to of the 
rheological flow curves measured in the presence of wall slip illustrated in Fig. 3(b) which consequently translates 
into a weak dependence of the hysteresis of the measured mean flow speed illustrated on to (see Figj7p. 

The central conclusion of the investigation of the wall slip behaviour is that, in the case of a unsteady low Reynolds 
number flow of a Carbopol® solution, the slip behaviour can not be decoupled from solid-fluid transition. Future 
theoretical developments should account for this experimental fact. 



4. Conclusions, outlook 

An experimental investigation of the slow yielding of a Carbopol® gel in both a unsteady rheometric flow between 
parallel and in a inertia-free unsteady pipe flow is presented. 



The first part of our study presented in Sec. 3.1 dealt with a detailed experimental investigation of the solid-fluid 
transition of a Carbopol® gel observed during unsteady flow ramps performed in a parallel plate geometry mounted 
on the rotational rheometer. 

The comparative investigations of the solid-fluid transition in the parallel plate geometry with and without slip 
performed for a wide range of characteristic forcing times to lead to the following conclusions: 

1. In the presence of wall slip, the solid-fluid transition occurs earlier, i. e. the apparent yield stress measured in 
the presence of slip t" is nearly an order of magnitude smaller than the yield stress measured in the absence of 



wall slip, T v (see Fig. 3(a)) 



In both the slip and the no slip case the transition from a solid behaviour to a viscous one (fully yielded) is not 
direct but mediated by an intermediate solid-fluid coexistence regime. The cusp visible within the solid-fluid 
coexistence regime on the decreasing branch of the stress ramp is associated to a elastic -recoil effect manifested 
by a flow reversal. 

In both the slip and the no slip cases the solid-fluid transition is not reversible upon increasing/decreasing of the 



applied stresses and a hysteresis is clearly visible (see Fig. 3(a)). 

4. The presence of wall slip nearly suppresses the scaling of the deformation power deficit P with the rate at which 
the material is forced, 1 /to- 

5. The irreversible solid-fluid transition and the wall slip are observed within the same range of applied stresses 
indicating that the two phenomena are coupled. 



The second part of our study presented in Sec. 3.2 concerns with a systematic experimental investigation of low 
Reynolds number unsteady pipe flows of a viscoplastic fluid (Carbopol®-980). By means of the DPIV technique a 
full characterisation of the flow fields during controlled stepped ramps of the driving pressure (Fig. || (b)) is made. 
Depending on the value of the applied pressure several flow regimes are observed. 

Corresponding to the lowest values of the driving pressure Ap a solid like flow regime is observed. This regime 
is characterised by a an elastic solid plug which never reaches a steady flow state. The largest level of fluctuations of 
the solid plug is observed near the centre line of the flow channel which we interpret as the signature of an unsteady 
solid (elastic) plug. 

For the largest values of the driving pressure Ap a fluid regime is observed. Within this regime the velocity 
reaches a steady state (see Fig. W) and the flow is re versi ble upon increasing/decreasing pressures, Figs, pi pi As 



in the case of the rheometric flow discussed in Sec. 3.1, the transition between the solid-like regime to the fully 



yielded regime is not direct but mediated by an intermediate solid-fluid coexistence regime. For each value of the 
characteristic forcing time to we have explored, this intermediate phase coexistence regime manifests itself by a 
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clear hysteresis (observed upon increasing/decreasing the driving pressure drop) which is qualitatively similar to 
the hysteresis observed in the rheological flow curves. The experimentally observed hysteresis of the flow curves 
is associated with a deformation power deficit and measurements of this quantity for various values of fo reveal 
no sensitive dependence on the frequency of forcing, A °c r ~° 1 (see Fig. Q). This fact is fully consistent with the 
measurements of the deformation power deficit performed in a rheometric (plate-plate) setup in the presence of slip, 
p oc f -° 03 , Fig. p(b). 



As a conclusion for this part, the main features of the solid-fluid transition observed in a rheometric flow in the 
presence of wall slip are equally observed in a unsteady pipe flow. 

In spite of its ubiquitous nature and tremendous importance, the phenomenon of wall slip in pasty materials and 
dispersions is still under-looked by many researchers, Buscall (2010). As Carbopol® gels are concerned, a detail 
characterisation of their slip properties in slow pipe flows was, to our best knowledge, insufficiently documented. The 
time and space resolved measurements of the flow fields based on the DPIV technique allowed us to attempt to fill 
this gap and provide a detailed experimental characterisation of the wall slip phenomenon. Systematic measurements 
of the slip velocity during controlled pressure ramps with various degrees of steadiness (i.e. for several values of 
the characterising forcing time fo) reveal a clear coupling between the irreversible solid-fluid transition and the slip 
phenomenon. A linear scaling of the slip velocity with the wall velocity gradients is found, Fig. ^. This scaling law is 
universal in the sense that it does not depend on the flow steadiness and invariantly holds for both the increasing and 
the decreasing branches of the controlled pressure ramp. 

A clear picture of the coupling between the irreversible solid-fluid transition and the slip phenomenon can be 
obtained by plotting the slip velocity against the wall shear stress, Fig. [h]. It is only within the fully yielded regime 
(A/? > Ap y ) that a Navier-type nearly linear scaling of the slip velocity can be observed. The measurement of the slip 
coefficient allowed one to estimate the thickness of the lubrication layer responsible for the slip behaviour. The result 
of this estimate seems to be consistent with previous experimental findings. It is equally interesting to note that this 
scaling of the slip velocity is practically independent on the flow steadiness (the characteristic forcing time fo) and 
whether the driving pressure is increasing or decreasing. 

Beyond the main findings briefly highlighted above, the present study also brings a modest contribution from 
a technical point of view: our approach of correlating DPIV measured flow fields with rheometric flow curve data 
allows a local measurement of the stresses which, in other words, "converts" the pipe flow into a table top open flow 
rheometer. To our best knowledge, this experimental approach is rather new (particularly in the context of viscoplastic 
flows) and deserves being pursued. 

We believe that our experimental findings could trigger several interesting developments at both experimental and 
theoretical level. A micro flow investigation of the wall slip of Carbopol® gels during unsteady flows over smooth 
solid surfaces able to simultaneously monitor the flow fields and the micro-gel structure near the solid surface would 
bring in our opinion a ultimate and direct clarification on the coupling between yielding and slip as well as on the 
microscopic scale mechanism of slip. Though we are aware that such experiments represent a daunting task, we note 
that several steps have already been t aken in this di r ection by combining microfluidics w ith various micro-rheology 
app roaches (diffusion of n ano tracers, Oppong et al. ( 2006 ); Dppong and de Bruyn ( 2007 ), fluorescent labelling of the 
gel, iGutowski et at] ( |2012| )). 

The success of such a task is strongly conditioned by a better understanding of the physicochemical properties 
of the Carbopol® gels which would allow developing versatile protocols of fluorescence labelling of the micro-gel 
particles which would allow the direct visualisation of the gel system under shear. 

From a theoretical perspective, our work clearly demonstrates the need of developing new frameworks to under- 
stand th e unsteady yielding a n d slip in Carbopol® gels. A modest progress in this direction has been already reported 
in Ref. Putz and Burghelea ( |2009| ) where a simple model that can accurately describe the unsteady yielding in a 
rheometric setup in the absence of wall slip was proposed and validated agains experimental data. 

In spite its ability to describe a large amount of rheological data in both she ar and oscillation and with in a wide 
range of polym er concentrations and at var ious operating temperatures(see Refs. Putz and Burghelea ( 2009| ); Moyers- 
Gonzalez et al. (201 1); Weber et al. (2012)), this model is minimalistic and can be considered for now as "poor man's 
model": it is a scalar model, it has not been derived from first principles, it lacks a thermodynamical validation and, 
in its present form, it can not account for wall slip. 

Future theoretical work is needed to correct these deficiencies and to ultimately attempt to model realistic vis- 
coplastic flows of industrial relevance such as flows around solid objects moving through a gel Putz et al. ( |2008| ), 
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unsteady start-up pipe flows and thermo convective flows. 
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